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We study the long-time aging dynamics of spin-glass models with two-spin interactions by perform- 
ing a Renormalization Group transformation on the time variable in the non-equilibrium dynamical 
generating functional. We obtain the RG equations and find that the flow converges to an exact 
fixed point. We show that this fixed point is invariant under reparametrizations of the time variable. 
This continuous symmetry is broken, as evidenced by the fact that the observed correlations and 
responses are not invariant under it. We argue that this gives rise to the presence of Goldstone 
modes, and that those Goldstone modes shape the behavior of fluctuations in the nonequilibrium 
dynamics. 
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I. INTRODUCTION 



Glassy materials are characterized by very slow dynamics, associated with a dramatic slowdown of molecular 
relaxation in structural glasses, and with a dramatic slowdown of spin relaxation in spin glasses. This slowdown of 
the dynamics has been captured in great part by the results obtained by dynamical mean field theories. In the case 
of supercooled liquids, the mean-field mode-coupling approach^ has been successful in describing some of the features 
of the relaxation. In the case of spin glasses, a dynamical mean field theory based on a Langevin dynamics for the 
spins, examined within a functional integral formulation of the Martin-Siggia-Rose approach 2 ^^! 5 -, has been used to 
study the long-time relaxation. The dynamical mean field theory of spin glasses has successfully captured 3 -! 4 ^ 5 - some 
unusual properties of the spin dynamics, associated with the lack of equilibration, including the presence of physical 
*^ ■ aging and the breakdown of the equilibrium fluctuation-dissipation relations. 

However, mean field theories do not allow direct access to a description of the fluctuations in the dynamics. It 
turns out that fluctuations in the dynamics of glassy systems can in fact be rather strong, as it has been underlined 
i ^ i ' by the discovery of dynamical heterogeneities^^. Dynamical heterogeneities are nanometer-scale regions of molecules 
rearranging cooperatively at very different rates compared to the bulk. Recent studies of material systems near 
their glass transitions have uncovered substantial experimental^ 9 - ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 and simulationa l 18 ' 19 ' 20 ' 21 ' 22 ev- 
idence for their presence. Various attempts at theoretically addressing these strong fluctuations have been made, 
involving, among others, the ideas of dynamic facilitatio n 23 ' 24 ' 25 ' 26 ' 27 , the presence of a "random first order phase 
transition" 28 : 29 ' 30 , or the use of diagrammatic methods to carefully re-analyze and extend mode-coupling theoi y 31 ' 32 . 
. However, a detailed theory that explains the dynamical heterogeneities remains elusive^ 3 -. 

Recently, a theoretical framework for the study of fluctuations in the non-equilibrium dynamics of glassy systems 
has been propose d 34 ' 35 : 36 ' 37 , which is based on the presence of a Goldstone mode associated with a symmetry under 
continuous reparametrizations of the time variable. It was argued there that the presence of this symmetry could 
provide an explanation for many of the dynamical heterogeneity effects observed in various glassy systems. In Rcf. 34 , 
L* ■ a sketch of a proof for the presence of this symmetry was presented. Earlier work had uncovered the presence 
of a restricted version of this symmetry, for the mean-field dynamical equations of some infinite-range spin glass 
models 3 -! 4 .^. 

In the present work, we present a detailed proof of the presence of this symmetry under continuous reparametriza- 
tions of the time variable, for the long time dynamics of a generic spin glass model with two-spin interactions. The 
proof is based on using the Renormalization Group to extract the long time behavior of the theory. It is somewhat 
unusual in the sense that we coarse grain time differences and not positions. In other words, the degrees of freedom 
that are "integrated over" are the ones associated with the "fast" dynamics, where by "fast" we mean fast in time, 
and not necessarily in space. 

Although involved in some of its details, our procedure is conceptually simple. We consider a model for a set of 
soft spins on a lattice, which contains only two-spin interactions, with a zero-mean uncorrelated Gaussian distribution 
for the spin couplings. We assume a Langevin-type dynamics for the spins, with a noise term whose amplitude is 
controlled by the temperature of the environment. We use the functional integral formulation of the Martin-Siggia- 
Rose approach to describe the Langevin dynamics. We set up the calculation by writing the generating functional for 
the spin correlations and responses, and find that this generating functional can be written in terms of a functional 
integral over an auxiliary field that depends on two times. We set up the Renormalization Group procedure by defining 
a cutoff To for the time differences. We increase the cutoff slightly, and integrate over all values of the auxiliary field 
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that correspond to time differences smalier than this slightly increased cutoff. This integral is actually a gaussian 
integral that can be performed exactly. After integrating over the "fast" variables, we rescale all times in such a way 
that the cutoff goes back to its original value To- We find that the RG flow converges to a fixed point, which defines the 
fixed point generating functional. Finally, we consider a smooth and monotonously increasing but otherwise arbitrary 
reparametrization of the time variable t — > s(t), which induces a transformation of the sources for the generating 
functional. We compute the value of the fixed point generating functional for those transformed values of the sources, 
and show that it is the same as for the original values of the sources. In other words, the reparametrization of the 
time variable leaves the fixed point generating functional invariant. 

The rest of the paper is organized as follows: in Sec. [II] we introduce and briefly discuss the spin model, the Martin- 
Siggia-Rose formalism for the Langevin spin dynamics and the assumptions about the nature of the random couplings; 
in Sec. IIIII we obtain an explicit form for the disorder-averaged Martin-Siggia-Rose generating functional, which 
contains the above mentioned auxiliary fields that play a central role in the formulation of the Rcnormalization Group; 
in Sec. IIVI we introduce our renormalization group procedure, associated with coarse graining the time differences, 
derive the flow equations for the parameters of the action, and find the fixed point to which the RG flows; in Sec. |V] 
we derive the central result of this work, i.e., we show that the fixed point generating functional is invariant under 
reparametrizations of the time variable in the sources; and in Sec. IVII we discuss the physical consequences expected 
from the presence of this symmetry, which have already been observed in numerical simulations of spin glasses and 
structural glasses, and can also be tested for in confocal microscopy experiments in colloidal glasses. Finally, in 
Sec. IVIII we summarize our results. 



II. MODEL 



We consider a spin-glass hamiltonian containing only two-spin interactions: 

H a = - Jrr'^Ar' +J2 W ^)- (1) 
TT t r 

Here the indexes r, r' label the N possible positions in the (discrete) lattice, the <f> r are soft spin variables, the J rr ' 
are the spin coupling constants (satisfying J rr > — J r > r and J rr ' = for r = r'), and the one-spin potential W(4>) is 
chosen to control the magnitude of the spin variables. We assume that the potential W((f>) is real, even and analytic 
at <j> = 0, i.e. 

oo 

W^)=5> p 4> 2p , (2) 

p=0 

with w p — w* Mp. For example, for the potential W{cf)) = j(l — (f> 2 ) 2 , the coefficients are wo = j, w\ = — ^, w 2 = j, 
and w p — Vp > 2. 

The Langevin equation for the spin variables for a given realization £ r (t) of the noise reads: 

We assume, as usual, that the noise is Gaussian distributed and uncorrelated, with a variance that defines the 
temperature T of the heat reservoir: 

(&(ti)&'(t 2 )) = 2T8 r ,r>5{ti-t 2 ), (4) 

where the angle brackets (• ■ ■ } indicate an average over the noise distribution. 
We compute the derivatives 



d(f> r 

where we have used that for all r, J rr = 0. 



a tt 

° = Y J Jrr^r' + W'{4> r ), (5) 
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Then the Martin-Siggia-Rose generating functional^ 9 -, averaged over the realizations of the noise, and incorporating 
the sources, reads 



(Z[i,h]) 



V<j)V4>V$ exp + Y J f f dt (tr(t)(f>r(t) + ih r (t)4> r (t) 



(6) 



where the notation ((•••)) indicates the average over the realizations of the noise. We are considering the time 
evolution between times to and tt of the spins 4> r (t), with initial conditions given by the ip r , i.e. Vr : <p r (to) = (f r , and 
the action is given (in general) by 



a 



In our case, by using Eq. ([5]) we obtain 



r J t 



•it) , 8H 



Of 



iT<f> r (t) 



(7) 



(8) 



We assume that the disorder is given by an uncorrelated, zero-mean, gaussian distribution for the couplings, i.e. 



(9) 



Here the connectivity matrix 2K rr i = J^ r , defines the variances of the random couplings, with the notation (• • • ) 
denoting an average over the disorder. The connectivity matrix K rr i encodes the properties of the model. For 
example, in the c ase of the Edwards- Anderson model, K rr > = K > for r, r' nearest neighbors and is zero otherwise. 

The functional (Z[t, h]) allows the direct computation of measurable quantities: expectation values, correlations 
and responses. The expectation values and the p-point correlation functions of the field <p r (t) are calculated by taking 
derivatives of (Z[£, h]) with respect to the source £ coupled to 



-(*)> = 



S(Z[£,h]) 



S£ r {t) 



'=0,h=0 



C p (ri,ti;r 2 ,t 2 ; ■ ■ ■ ;r p ,t p ) 



n{tl)(pr 2 (t2) ■ ■ ■ 4> rp {tp)) 



S^(Z[£,h]) 



5£ ri (ty r2 (t2)---£ rp (t p )- 



(10) 



(11) 



Here we have used the fact that the generating functional (Z[£,h\) reduces to unity for zero sources, i.e. it satisfies 
the condition (Z[£ = 0,h = 0]) = 1. 

The effect of (possibly time dependent) external fields: H = H — J2 r h r {t)<j) r {t) can also be probed by computing 
response functions: 



R(r,t\r',t') 



5(Mt)) 
5h r >[t>) 



6{Z[£,h}) 



e=o,h=o 



5£ r (t)5h r ,{t') 

i((f>r(t)$r>(t')}, 



X (r,t\r',t') = / dt"R(r,t\r',t"). 



(12) 
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Here R(r,t\r' ,t') represents the response to an external field only present at time t', i.e. a "delta function in time", 
and the integrated response x(( r > A r ' 1 O corresponds to a "step held", i.e. an external held the is "turned on" at 
time t' and "stays on" until the time t when the spin is measured. An important property of response functions 
is that, by causality, the response R(r,t\r' ,t') is zero for t' > t. Expectation values, correlations and responses for 
one disorder realization can also be computed by formulas that differ from Equations (jlll) and (| 1 2[) only in that all 
disorder averaging is removed. 

III. DISORDER AVERAGED GENERATING FUNCTIONAL 

Once the distribution of the couplings is defined, we can average the disorder-dependent exponential in the action: 



Z (/>,< 



exp —i / dt J rr >(j> r (t)(f> r > (t) 

\ r,r' 

exp(-J r 2 r ,/4^, 



rr' I 

X 



= n / dJr 

r<r' 

exp (^J rr i |— i J dt (j> r (t)(f> r i(t) + (j> r i(t)4> r (t) 

= exp < - V] ^TT~ f dtxdte (<f>r(h)<f>r' (h)4>r(t2)<t>r' {h) + 4>r>(tl)4>r(tl)4>r(t2)4>r>{t2) 
[ r,r< 1 Jt V 

+^r(ti)^ r /(ti)^(t 2 )^r(t2) + ^r'(*l)^r(*l)^r'(*2)^r(<2))}. (13) 

We now dehne the notations = <j) r (t), 4>\{t) = (f) r (t), = 1, 1 = 0, which allow us to write: 

Zo[$,<f>] = e X pl~J2 R rr' dhdt 2 J2 MtlWW&itlWto)) ■ ( 14 ) 

Y r,r' ^ to a : c=0 J 

Here we can introduce auxiliary two-time fields Qr C {h,h) by performing a Hubbard-Stratonovich transformation: 
z [$,<j>] = (VQ exp [ ~^M rr , I dhdh ^ QTih,h)QlKh,h) 

J \ r,r' ^ to a,c=Q 

+tJ2 f h dhdt 2 ]T Q™(h,t 2 Wr(h)<t> C r (h)) > (15) 



where M rr , is the matrix inverse of K„, and J VQ = M{M) J H r a c Ut u to dQ^ c (h,h). Here N{M) 

1/2 

(det ((27r)~ 1 Af)) is an M-dependent normalization factor. 
We're now in a position to write down the disorder-averaged generating functional for the problem: 



Z[£,h] ee (Z[£,h]) = / VQ eiq>(-S K [Q]-S nl [Q,£,h]), (16) 
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where 



S K [Q] = \ Y. M " J f 1 dhdt 2 J2 Qr C (ti,t 2 )Q^(h,t2) (17) 

r,r' •' t ° a,c=0 

S„i[Q,£,h] = -In jvcjPV^Vy exp {iS HS [Q, + iS spi n[<t>° , t 1 } + iS BC [^, <p] 

+ E I 1 dt Ur(t)(f) r {t) + lh r (t)j> r (t)] 1 (18) 
r Jto ) 

ShsIQ,^^ 1 } = P dhdh E Q a r C {tiM)€ihW r (t 2 ) (19) 

r J *<> a,c=0 

St&iM = -e jytm^^+w'i^-iTm) 



(20) 



S BC [<i>\0\ = X)^K(*o)-Vr}. (21) 

r 

To simplify the algebra, we take from now on the integration limits as to — and tf — > oo. By combining Eqs. ([2]) 
and (120|) . we write S^pm in a way that will allow the RG equations to be put in a simple form: 

SsphJ^ 1 ! = -V f°°dt ( ^Sft)S+ v ^rjm^+T^U^micnr 1 - 1 



■°,^] = -E /"* ff^w^lr + E 7i?cw^(i)+E7 (2rt ^w(^w) 2 

/ dt x dt 2 y^K rr , E fl (4) (*i-*2)^(*i)^(t 2 )^(*i)^'(*2). (22) 



a,c— 

IV. RENORMALIZATION GROUP 

We want to introduce an RG transformation on the time variables. Since the construction of the RG transformation 
is a bit unusual, we will explain it in detail. We introduce a short-time cutoff tq = l/fio for the time difference 
t\ — t 2 . This only affects the terms in the action containing an integration over two time variables, namely Sk[Q], 

spin [ 



Shs[Q, 0°) ty 1 ] an d Sspin [4> , 4> ]■ The first two terms take the following form as a starting point for the RG 



1 r 1 

S K [Q] = oE M -' / dt ^ E QT{tiM)Qf{txM), (23) 

r,r' o<t 1 ,t 2 <oo a,c— 

ShsIQ,* ,* 1 ] = E / E Qr%h,h)<t> a r{tl)<t>r(h)- (24) 

/ T 0< 1*1 -'2 I 

r o<t 1( t2<oo a,c— 

These terms differ from Eqs. (fT7| and (fl"9|) by the removal of the contributions corresponding to \t\ — t 2 \ < t . 
There are two possible natural assumptions about how this cutoff is implemented: either we assume (i) that the 
contributions for those time pairs is directly removed from Sk[Q] and Shs[Q, <fi°, 1 ] without any effects on other 
terms in the action, or (ii) that the Hubbard-Stratonovich transformation performed to obtain Eq. (|15p is undone for 
time pairs \t± —t 2 \ < tq. These two alternative assumptions lead to slightly different starting points for the RG flow, 
but in the end the flow converges to exactly the same fixed point in both cases. This is reassuring, in the sense that 
we expect the properties of the long time dynamics not to depend on the cutoff procedure. The initial coefficients for 
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Sspin^ , 4* 1 } in Eq. are given by 



r = 1 

(2) -rr 

(2) (2) r, 

7oi = Tio = 2 ™i 

v( 2 P) 



= 2pw p Vp>2 (25) 

(4) , t ) — I ® cuton " P roce dure (i) , . 

^ 1 2 > ~~ \C\t 1 -t 2 \<T Q for cutoff procedure (ii) . ^ ' 

Here the characteristic function C-p is defined to be 1 if V is true and if V is false. 

We now perform an RG transformation on the time variables. We separate the two-time fields Q into fast modes 
Q > and slow modes Q<: 

n ac j. \ _ / Qr c (ti,t2) for T < \h - t 2 \ < br ( 
^A^) = \ for 6r <|ii-i 2 |, { ' 

n ac ( . . s = J for r < |*i - 1 2 | < &t , , 

Li< ' Atllt2 > ~ \Q a r c {hM) for 6r <|ii-i 2 |, [ ' 

with 6 > 1. Clearly, we have 

Q^(t u t 2 ) = Q^ r (t 1 ,t 2 ) + Q a < c < r(h,t2), (29) 
and by inspecting Eq. (|23p we find that 

= 5k [Q> + Q<] = 5k [Q>] + S*[Q<]. (30) 
As our next step, we integrate over the fast variables Q > to obtain 

Sn/ b [Q<] = -In/ 2?Q> exp{-5n[Q>+Q<]} 

1 /" 1 

= F n + -J2 M rr' / dtidta V Q a <M 1 ,t2)Q< r ,(t 1 ,t 2 ) 

Z * — ' 7 &to<I*i— *al ^ — ' 

r,r' o<t l7 t 2 <oo a,c— 

-In / VjPVtfVQ \ f PQ>exp ( -\ Y^M rr , [ dhdt 2 V Q§.%(ti,t 2 )Q§fv(ti,*a) 

/ / I Z ' / to < l*i -*2 1<<"- ' ' 

17 l \ r ' r ' 0<t llt2 <oo a,c=0 

+*E / o , ^1*2 ]T Q£(*i,k)#(*i)#(*a) )) 

* ' / ^0<l t l- t 2l< b ' r * ' / 

r o<t 1 ,t 2 <oc a.c— / ) 

xexpfzV / dt 1( ft a V Q^ ^ ,(^l,^ 2 )C(^l)^(^ 2 ) + iS , S p in [ ( /) o ,0 1 ]+^5 B c[0 1 ,^ ) (31) 

\ ' ' / ("T <|t 1 -i 2 l * ' / 

\ r 0<ti,t 2 <oc a,c— / 

The factor {J'DQy cxp{- • • }}, which contains the integration over the fast modes Q>, is actually a gaussian 
integral, which evaluates to 

(det{(2 7 r)- 1 ^}) V(T0 ^ )/2 

xex P(-^E X -'/ dhdhY] ^(*i)^(* 2 )^(ti)^(f2) | , (32) 

\ 2^^ /to<I*1-*2|«'t / 

\ r,r 0<ti,t2<cc ac— / 



where V(ro, 6) is proportional to the volume of the two-dimensional (time) region where the condition tq < \ti — 1 2 | < 
&ro holds. In this expression, the determinant prefactor contributes to the renormalization of the constant term Fq, 
and the argument of the exponential contributes to the renormalization of the function g( 4 >(ti — t 2 ). 
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We now not only perform the rescaling of the fields Q< c r (ii, tz) and the time variable, as it would normally be done 

for an RG procedure, but we also simultaneously rescale the <f>®(t) fields and the sources {£ r (t), h r (t)}, even though 
those quantities were not subject to the integration of "fast modes" : 

QZr(KM' 2 ) = ^ Q'riAA) (33) 

bt' = t (34) 

#(6f) = b^&it') (35) 

£ r {bt') = b Xe e' r (t') (36) 

h r {bt') = b Xh h' r {t'). (37) 



We then get 



s'AQ'} = (& 2+A « +A «) l^M^ L, t ,_ t ,, dAdt' 2 Y, Q' r a \AA)Q? c {AA) (38) 

r,r' o<t^,tj<oo a,c=0 

S' HS [Q'J ^] = (b^+^+^Y L lt ,_ t , } dt[dt> 2 Qr{AA)l a {A)l c {Q- (39) 

r 0<t[,t' 2 <oa a : c=0 

Since the terms Sk and Shs together represent the 4-spin interaction that makes the model glassy, we demand 
that they both should be marginal under the RG. This leads to the conditions 

= 2 + (40) 
= 2 + Ai 2 c >+A( 1 )+A( 1 >. (41) 

The second condition can only be satisfied if is of the form A^ = /(a) + /(c), where f(a) = — 1 — A« . Inserting 
this form into the first condition, it yields 1 + f(a) + f(a) = 0. At this point we still have freedom to pick among 
infinitely many possible solutions to this equation, each one of them defining a different RG transformation. We 
decide to choose the assignment f(a) = —a, which leads to 

A£) = -a-c (42) 
Ai 1} ee a - 1 = -a. (43) 

The choice of this particular solution is natural if we consider a Reparametrization Group (RpG) transformation^, 
associated with a reparametrization s(t) of the time variables, 

QTit'xA) = (Jf)" (j^rj : Q a r c (s(t\), s (t> 2 )), (44) 

where a,c £ {0, 1}. For the special case of a rescaling of times, s(t') = bt', Eq. (|44]) reduces to 

Q a r c (bt' 1 ,bt' 2 ) = b- a - c Q a r c (t , 1 ,t , 2 ), (45) 
which is completely analogous to Eq. ()33j) in the case Xac — & c. 



For the source term, we demand that it should be marginal under the RG, and obtain the rescaling exponents: 

A, = -1 - A^ = -1, (46) 
A,, = -1 - A^ 1} = 0. (47) 

It can be checked that, besides Sk, Shs and the source term, the boundary condition term Sbc is also marginal 
under the RG. 

We now consider the effect of the RG transformation on the terms contained in S sp i n [(f> , (j) 1 ]. The time derivative 
term is not affected by the integration over fast modes, and the rescaling of times and fields introduces the following 
rescaling: 

1 1_ _ _ J_ 



8 



If we now write 



we get the RG equation 



b = e sl , (49) 



Similarly we obtain 



(2) 

ac 



dl 

d7 (2p) 



dl 



dr , x 



= ( + c -l) 7 W, (51) 
= 0. (52) 



Finally, from Eq. (|32p . we find that the only term to receive a contribution from the integration over the fast degrees 
of freedom is the g^{t\ — ta) term: 



(h - ta) -> .9 (4) (^'i - 4) = 6 2 - a - £ - a - c x (>)(&*; - M' 2 ) + C ro <| fct; _ ht , , <bro ) , (53) 



and we observe that for this term the rescaling prefactor evaluates to unity: b 2 a c a = 1. 

By examining the RG flow of Eqs. (|50"|) . (j5Tj) . (j52|) , and ([53]) . we hnd the following fixed point values: 



T* = 0,oo 

*(2) n 

7oo = 0, oo 

= any number 
7io^ = any number 

*(2) n 
7ll = 0,00 

^,*(2P) _ nUm b er \fp > 2 

5* (4) (^i-^2) = C| tl _ t2 | <ro . (54) 
Since the RG flows of all parameters are uncoupled, the solutions above can be chosen independently for each param- 

(2) (2) 

eter. The stability analysis around the fixed points shows that perturbations of T, 7 q , and 7ll are relevant near 
r* = 0, 7qq 2 '' = oo and 7^ 2 '' = respectively; and are irrelevant near T* = oo, Jqq 2 ^ = and 7*} = oo respectively. 
It also shows that perturbations of 7qi\ 7io\ and ^ 2p ^ are marginal around any of their fixed points. Perturbations 
of g( 4 '(ti — £2) are always irrelevant. Therefore, for the set of initial conditions given by Eq. (f2T)|) , and for almost any 
other set of initial conditions, the RG flows for T, 7qq\ and g^(ti — £2) converge to their stable fixed points. However, 
the parameter 7 X1 has a starting value which is exactly at the unstable fixed point = 0, and stays there through 

the RG flow. Additionally, the parameters 7qi\ and j@ p ") do not flow at all, and stay at their initial values. In 
summary, the parameters of S sp - m [cjP , (j) 1 ) flow to the fixed point values: 



J1* 


= OO 


*(2) 

7oo 


= 


*(2) 

7oi 


= 2w 1 


*(2) 

7io 


= 2wi 


*(2) 

7n 


= 


7 *(2p) 


= 2pw p Vp > 2 


g* (i) {t 1 -t 2 ) 


= C\tl-t 2 \<T - 



(55) 

As anticipated above, this result is the same for cutoff procedures (i) and (ii) (and in fact for any other possible initial 
value of gW(ti -t 2 )). 

The fact that T flows to infinity indicates that the derivative term does not appear in the fixed point action. 
However, the states of the system at different times are still coupled by three other terms: Sk[Q], Shs[Q, 4>° ■> and 
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the term proportional to g^\t\ — t 2 ) in S'spinfc/) , We interpret this to indicate that, while the time derivative 
terms may be important for the short time dynamics, when the short time dynamics is "integrated over" and only 
the long time dynamics remains, the coupling between different times is provided only by the terms associated to the 
spin glass interactions. This is reminiscent of earlier mean-field calculations of the aging dynamics of spin glasses, in 
which the time derivative terms are negligible at long times, and the coupling between different times is also provided 
only by the spin-glass interaction terms^. In that context, the time derivative terms break the mean-field version of 
time reparametrization invariance, and therefore this invariance is only valid for very long times, when time derivative 
terms are negligible^^^. 

By combining all the results for the RG flow for the various terms in the action, we find that the action converges 
to the fixed point: 

1 r 1 

5 fp [Q] = Fn + -J2M rr , / dhdt 2 V Q a r c {tiM)Qf{tiM 

Z * ' / T)<|tl -*2 I £ * 

r,r' o<t 1 ,t 2 <oc a.c—0 

-In [v<tPV<l> l Vy exp/iV / dhdt 2 V Q? c (ti,t 2 )0?(ti)^r(*a) 

J J T <\t 1 -t 2 \ 

\ r o<t l7 * 2 <oo a,c=0 
poo I oo \ 

- j E j o dt (-2^i^W^(t)-E 2 f^^(*)(^(*)) 2p " 1 J 
i r 1 

-iY,Ktv> / dhdt 2 V #(ti)#(t 2 )#(*i)#(*a) 

I ^— ' J |ti-*2l<*no 

r,r o<*i,t2<oo a,c=0 

+»53v?r{^(0)-^r} + 53/ * (4(t)^(*)+tM*)0?(*)) [• ( 56 ) 



In this form the fixed point action no longer contains the auxiliary fields Qr C {ti,t 2 ) f° r times ti,t 2 such that 
|ti — ^2 1 < tq. We now re-introduce those auxiliary fields through the same Hubbard-Stratonovich transformation 
that was used to obtain Eq. ([14]) . and obtain: 



S fp [Q] = Fn + i^)^ // dt x dt 2 Qr C (ti,h)Q*?(h,t 2 

r.r' a,c=0 



In / V()PV4> l V$ expliV f f dt x dt 2 V Q^ii, ta)#(*i)#(*a) 

/.OO / \ 

-*E y o d * ^-2^i^(t)^(t)-^2p Wp ^(t)(4(t)) 2 f- i J 

+i^2tPr{4>l(P)-<Pr} + ^J di (M*)#(*)+*M*)#(<)) [• (57) 



V. REPARAMETRIZATION SYMMETRY 



We are now finally ready to evaluate the effect of a reparametrization i — > s(t) of the time variable on the fixed-point 
generating functional Zf p [£, h]. We consider any smooth monotonous increasing function s(t) satisfying the boundary 
conditions s(0) = and s(oo) = oo, which induces the following transformation of the sources: 

£•(*) = (SV(«(*)) (58) 



ft P (t) = M«(t)), ( 59 ) 
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and compute the fixed point disorder averaged generating functional, evaluated at the transformed sources: 
Z ip [{£ r (t),h r (t)}} = 



J VQ CX P j -F U -~J2 M rr' J J dhdt 2 Qr C {h,t 2 )Q^{tlM) 
\ r,r' a.c—0 

+ ln f VjPv^VQ exp/zj] 11 dt x dt 2 V Q a r c {t x MW T {txW r {t-i) 

J { r J J a ,c=0 

poo / oo \ 

/ dt \ -2^i^(«)^(*)-E 2 f w v tfitM®) 3 *- 1 

r J \ p=2 ) 

+i^2<p r {lpl(Q) - Lp r } 
r 

+ J2 f dt (4(t^(t) + ih r (t)i>° r (t)) I j . 



(60) 



Here we have changed the name of the dummy variables from (p to ip and from Q to Q m the functional integral. 
We now perform the changes of variables ip$ it) = (§§) a <^(s(i)) and Q^ c (h,t 2 ) = (^7)° (j^Y Q£ c (s(*i), s(t 2 )), i.e. 
the change of variables associated with the RpG transformation of Eq. (|44[) . thus obtaining: 



z {p [{e r (t),h r (t)}] 

PQ «'xi» | -F n - i^M rr / / / dt x dt 2 



a,c— 



5s 



a+a 



<3s 
9^ 



+ ln / V(f>° Vc/) 1 Vip expii^ / / diidt 2 



5s 
9^ 



a+a 



ds_ 



xQrWi)^(i2))^(s(«i))^(s(t 2 )) 
4(8(t))4(a(t)) -p2P w P (^) OtfOK*))) 2 '" 1 ) 



+*X)^r{0j(s(O))-V>r} 
r 



2? I 



AO ,A1 



lnj 2 



Here the symbol J7i 



represents the Jacobian of the transformation from cp to ip, and the symbol J 2 



(61) 



2?Q 



represents the Jacobian of the transformation from Q to Q. Since both transformations are linear transformations, 
the Jacobians only depend on the reparametrization s(t), but they do not depend on the fields <p or Q, or the sources 
{£ r [t), h r (t)}. For this reason, we will denote them as Ji{s(t)} and J 2 {s(t)} respectively. Using the fact that a + d — 1 

, \ a+a / „ \ c+c 

and c + c = 1, we find that the factor I J ( j^- J is simply the Jacobian of the transformation from [tx , t 2 ) to 
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{t'lit'?) = (s(ti), s(t2)), and therefore we obtain: 
Z ip [{I r (t),h r (t)}} 

= JvQ exp j -F a -^ M ^ / dt t dt 2 E Wi.4)<9r"(*i 5 4) 

\ r,r' a,c— 

+ ln y VcfPVcfVv exp!iJ2 j ^ ' « £ Q? c (*i>*2) 



a,c— 



poo / oo > 

/ dt ' -2t 01 ^(o^(t , )-E 2 p^(0(^(*')) 2p " 1 

+i^9? r (<#(0) -yv} 



+ E / o d *' (4(t')^(^)+^(i')^(0)+lnJi{sW}| 
+ lnj 2 {s(t)}) 

= Z tp [{£ r (t), h r (t)}} x Ji{s(i)} x J 2 {s{t)} (62) 

Here we have used the boundary condition s(0) = 0. We now consider the special case of zero sources, i.e. t r (t) = 0, 
h r (t) = 0; in this case, the transformed sources are identical to the original ones, and we have the condition 

Z fp [{0,0}] = Z fp [{0,0}] x Ji{s(t)} x J 2 {s{t)}. (63) 

Since the generating functional is nonzero for zero sources (it is actually unity^ 9 -), we immediately conclude that, for 
any reparametrization s(t) the product of the Jacobians is unity: Ji{s{t)} x ^{s^)} = 1. Thus we obtain, for any 
reparametrization s{t), the identity: 

Z fp [{l r (t),h r (t)}} = Z lp [{£ r {t),h r (t)}], (64) 

i.e., we have shown that the fixed-point generating functional is invariant under time reparametrization transforma- 
tions. 



VI. PHYSICAL CONSEQUENCES OF THE TIME REPARAMETRIZATION SYMMETRY 

Since the renormalization group procedure described above involves integrating over all short time-scale fluctuations, 
the fixed point generating functional that we obtained controls the long-time dynamics of the model. The group of 
transformations associated with time reparametrizations is a continuous symmetry group for the fixed point generating 
functional. This symmetry is broken by the actual dynamical correlations and responses observed in the system. As 
an example, let us consider the space-averaged two-time correlation C(t,t w ) = jj^ri't'ri^^ritw)), where t w is 
normally referred to as the "waiting time" and t as the "final time" . If this correlation was actually invariant under 
time reparametrizations, we would have C(t,t w ) = C(s(t), s(t w )) for any arbitrary increasing function s(t) such that 
s(0) = and s(oo) = oo. The only possible way that this condition can be satisfied is if C(t,t w ) = Co (a constant). 
Since correlations in spin glasses actually do change with time, this implies that the reparametrization symmetry must 
be broken. 

We have therefore the presence of a broken continuous symmetry group. Since no long range interactions or gauge 
potentials are present, we should normally expect that a Goldstone theorem applies, giving rise to the presence of 
Goldstone modes (or soft modes) in the system^. For this reason, it has already been argued in Refsi 34 ' 35 i 36 ' 37 that 
Goldstone modes should be present in the non-equilibrium dynamics of spin glasses and possibly other glassy systems, 
and could in principle constitute the main source of fluctuations in the non-equilibrium dynamics of these systems. In 
other words, the presence of time reparametrization symmetry could account for a significant part of the dynamical 
heterogeneity effects observed in glassy systems. 

In general, Goldstone modes are obtained from a continuous symmetry transformation by making it smoothly space 
dependent. For example, if the symmetry corresponds to spin rotations by any angle, to obtain a Goldstone mode 
the angle is chosen to be smoothly space dependent. In the present case, the continuous symmetry corresponds to 
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reparametrizing the time variable t — > s(t). In the uniform case, this leads to the symmetry transformation C(t, t w ) — > 
C'(t) *tu) = C (s{t) , s{t w )) . The Goldstone modes are obtained by choosing the time reparametrization to be smoothly 
space dependent, i.e. t — > s r (t), and C r (t,t w ) = Co(s r (t), s r (t w )), where Co(t,t w ) is space independen t 35 ' 36 ! 37 . Since 
the reparametrization is non-uniform, it is no longer a symmetry transformation for the system, but if the space 
variation is slow enough, the change in the action with respect to the value for a uniform two-time field is small. A 
possible (very simplified) physical interpretation of these Goldstone modes, is that they are associated with "non- 
uniform slow relaxation" : if one considers different small regions in the system, for all regions the relaxation path is 
very nearly the same (as given by Co (t,t w )), but the rate at which each small region advances in its relaxation path 
can fluctuate from region to region. 

Testing for the presence of fluctuations associated with this reparametrization symmetry has been performed in 
numerical simulations of both spin glasses and structural glasses. In spin glasses, one prediction that can be tested 
in simulations refers to the values of coarse grained local correlations C r (t,t w ) = ^ X^gb s i(t)si(t w ) and integrated 

responses Xr{t,t w ) = f* dt' ^J2ieB r S sh*(t>) > wnere Br is a cubic coarse graining box containing n spins centered 
at the point r; in comparison to the global values C(t,t w ) and x(t,t w ) obtained by taking the averages over the 
whole sample. As explained in Refs i 35 ' 36 , the presence of Goldstone modes associated with time reparametrization 
symmetry would imply that the pairs (C r , Xr) should be concentrated predominantly along the parametric curve x(C). 
It turns out that this is exactly what is observed in the results of numerical simulations in the 3D Edwards- Anderson 
mode l 35 ' 36 . Another testable prediction is that, if the global correlation C(t,t w ) is only a function of the ratio t/t w , 
i.e. C(t,t w ) — C(t/t w ) , the probability distribution p(C r (t,t w )) for the values of the local coarse grained correlation 
C r (t, t w ) should collapse as a function of t w , as long as t/t w is held fixed. This has also been found to be the case in 
simulations in the 3D Edwards- Anderson mode l 35 ' 36 . In Refill a more detailed study of the shape of the probability 
distributions for both the Edwards- Anderson model and a kinetically constrained model of glassiness was performed, 
with results that were consistent with the predictions derived from the presence of Goldstone modes in the system. 

Another aspect of the results presented here that can be tested by comparison with numerical simulations in spin 
glasses is the fact that the symmetry is only exact for the fixed point generating functional, i.e. in the limit t — ► oo. 
For the case of an exact continuous symmetry, one should expect that the presence of a true Goldstone mode (with 
zero mass) gives rise to spatial correlations that decay as power laws at long distances. However, at any finite time 
the symmetry is broken by small corrections to the action, which we can think of as small symmetry-breaking fields 
that go to zero at t — > oo. As a consequence of the presence of these symmetry breaking fields, the Goldstone modes 
now acquire a small mass, which should vanish in the t — > oo limit. In Ref.— , the spatial correlation length £(t,t w ) for 
fluctuations of the quantity Q]}{t, t w ) was measured in large scale long time simulations in a 3D Edwards- Anderson 
model. For very large t, t w and t/t w the time dependence of £(t,t w ) was found to be consistent both with a form 
£(t,t w ) ~ ln(tt w ) or a form £(t,t w ) « (tt w ) a , with a w 0.04. Both forms extrapolate (albeit slowly) to infinity at 
infinite times. This is suggestive, and consistent with what is expected from the results of the present work, but the 
actual values of £(t, t w ) are too small to make any firm statements about the t — > oo limit. 

The present work only proves the presence of time reparametrization invariance in spin glasses. However, it is 
conceivable that the symmetry could extend to structural glasses, and there has already been some work in structural 
glasses which has found suggestive evidence for its presence. In the case of structural glasses, there is another quantity 
which plays the role of local coarse grained correlation. It is defined^ as C T (t, t w ) = N ^ B < Y^r-(t )es cos(q • (ij-(i) — 
Tj(t w ))). Here B r denotes a coarse graining box centered at the point r in the system, and the sums run over all of 
the N(B r ) particles present in B r at the waiting time t w . The value of q is usually chosen to correspond to the main 
peak in the structure factor S(q) of the system. Unlike in the 3D Edwards-Anderson model, in structural glasses 
the global correlation C(t,t w ) is not a function of the ratio t/t w . For this situation, the presence of the Goldstone 
mode associated with time reparametrization invariance implicates that the probability distribution p(C r (t,t w )) for 
the values of the local coarse grained correlation C r (t,t w ) should collapse as a function of t w , as long as the global 
correlation C(t, t w ) is held fixe d 36 ' 37 . This has been found to be the case, to a good approximation, in simulations in 
binary Lennard- Jones mixtures and binary Weeks-Chandler-Anderson mixtures^. 

Confocal microscopy experiments in colloidal glasse o 9 ' 10 ' 11 provide detailed data that include the positions of all 
colloidal particles in some sub-volume of the sample at different times in the evolution of the system. These data 
can be analyzed in completely analogous ways to those used to analyze data from simulations in structural glasses. 
It remains an open question whether or not such analysis would provide further evidence in favor of the presence of 
time reparametrization symmetry. 
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VII. SUMMARY 



In this work, we have presented a detailed proof of the presence of a symmetry under continuous reparametrizations 
of the time variable, for the long time dynamics of a generic spin glass model with two-spin interactions. No assump- 
tions were made about the range of the interactions, therefore the proof applies equally to short-range models, such as 
the Edwards- Anderson model, and to long-range models, such as the Sherrington-Kirkpatrick model. By performing 
a Renormalization Group procedure that exactly integrates over degrees of freedom associated with short time differ- 
ences, we have obtained the RG flow for the parameters in the action. We have found that the RG flow converges to 
a fixed-point generating functional, and we have explicitly written the form of this generating functional. Our main 
result is to have shown that the value of the fixed point generating functional is left invariant by a transformation of 
the sources induced by a monotonous increasing but otherwise arbitrary reparametrization of the time variable. 

The group of transformations associated with time reparametrizations is a continuous symmetry group for the fixed 
point generating functional. This symmetry is broken by the actual dynamical correlations and responses observed 
in the system. In a situation like this, one would normally expect the presence of Goldstone modes. Indeed, it 
has been argue d 34 i 35 ' 36 i 37 that Goldstone modes associated with time reparametrization invariance should dominate 
the fluctuations in the non-equilibrium dynamics of these systems. Positive evidence for this statement has been 
found in simulations of the aging dynamics of the 3D Edwards- Anderson mode l 35 i 36 ' 37 . Even simulations in systems 
without quenched disorder, such as kinctically constrained models of glassiness^ 7 - and models of structural glasses 22 
show evidence in favor of the presence of this symmetry. Additionally, experimental tests for the presence of this 
symmetry in colloidal glasses can be provided by confocal microscopy measurements. Having proved the presence 
of time reparametrization symmetry, the present work opens the door for a more detailed analytical study of the 
symmetry itself, of the Goldstone modes probably associated with its presence, and more generally of the fluctuations 
("dynamical heterogeneities") that are present in the slow dynamics of spin glasses and other glassy systems. 
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